A gluon condensate term in a heavy quark mass 
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We investigate a connection between a renormalon ambiguity of heavy quark mass 
and the gluon condensate contribution into the quark dispersion law related with a 
virtuality defining a displacement of the heavy quark from the perturbative mass- 
shell, which happens inside a hadron. 
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An Operator Product Expansion (OPE) is among the most powerful tools in the heavy 
quark physics. In this respect it is usually applied in the form of series in the inverse 
heavy quark mass, determining the characteristic energy scale, say, in sum rules or for 
decays etc. It is well recognized that the Wilson coefficients standing in front of quark- 
gluon operators can contain the uncertainty caused by the factorization of perturbative 
contribution and the nonperturbative matrix elements of composite operators. In this case 
the restriction on internal virtualities in Feynman diagrams has to be introduced to control 
the dependence on an "infrared" energy scale A. Usually, the gluon propagator is modified 
by replacement: Xjk^ V(^^ ~ ^g) ^^e cut off the gluon momenta is performed as 
A;^ > 2]. The calculation results depend on these parameters. Say, a peculiar behaviour 
at A^ — i> appears in physical quantities. For example, a perturbative correlator of two 
heavy quark currents acquires a power correction like A^/m^, where m is the heavy quark 
mass jj]. Physically, it means that the OPE can be valid if we sum the perturbative and 
nonperturbative parts with the vacuum expectation of gluon operator which has the same 
low energy scale dependence: the gluon condensate ~ A^. Then the A-dependent term can 
be adopted by an appropriate definition of OPE with the condensates. Another case takes 
place for the uncertainty in the heavy quark mass, where the perturbative calculation of 
self-energy with the gluon virtuality cut off leads to the linear term in A. However, there 
is no appropriate operator whose vacuum expectation is proportional to the first power 
of low energy scale 1|. It was shown that the mentioned uncertainty proportional to the 
powers of factorization scale A can be related with the perturbative summation of higher 
order diagrams, which in the limit of infinitely large number of fiavors has the divergency 
of series in /^oa^, where /3o denotes the first coefficient of Gell-Mann-Low function in QCD. 
The Borel transform of such series has some peculiar points, which provide the uncertainty 
in the inverse transformation. This uncertainty, related with the divergency of perturbative 
series is called the renormalon |^, since the physical contents of such fact is clarified by 
the representation, where the series are combined in the running of QCD coupling constant 
dependent of the gluon virtuality. The coupling has the singularity, which is the indication 
of confinement. In this way, the uncertainty in powers of Kqcd appears again. Modern 
studies on the renormalon applications can be found in . These facts imply that the OPE 
for fixed values of physical quantities (say, partial widths or coupling constants in the sum 
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rules) in terms of perturbative heavy quark mass results in the heavy quark mas, whose 
value extracted form the data, strongly depends on the order of calculation in a^-series 0: 
the mass value is significantly changed from order to order. 

Thus, the heavy quark quantities have the renormalon uncertainties connected to the 
infrared confinement in QCD. Some of them can be eaten by the appropriate definition of 
OPE with condensates. The heavy quark mass is of a special interest, since its infrared 
uncertainty cannot be straightforwardly adopted by the vacuum expectation of an operator 
with the dimension 1 in the energy scale. 

In present paper we evaluate the gluon condensate contribution to the dispersion law of 
heavy quark. We find that the corresponding operator is divided by the third power of quark 
virtuality, which results in the appropriate dimension of term in the heavy quark mass. We 
discuss how this fact can be used to cancel the infrared uncertainty of mass. 

We perform the calculation of diagram shown in Fig. 1 in the technique of fixed-point 
gauge p] with the NRQCD propagators of heavy quarks 0]. 
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FIG. 1: The diagram with the gluon condensate contribution to the two-point effective action of 
heavy quark. 



The covariant form of two-point heavy quark effective action hyTh^ can be represented 
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where v denotes the four-velocity of hadron containing the heavy quark. The validity of 
holds under the certain condition on the region of kinematical variables: the gluon 
condensate term in the dispersion law of quark is less than the leading contribution. 
In the rest frame of hadron v = (1,0) we have 
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where Ai? denotes a heavy quark virtuality inside the hadron. The perturbative mass-shell 
is defined by the following expression: 



AE = 0. 
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It is quite clear that the confined quark cannot reach the mass-shell and there is a minimal 
displacement from the surface of free quark motion, which is a nonperturbative quantity. 
So, we suppose that 

In what follows we apply the model with the quark dispersion law determined by the form 
dictated by the account of gluon condensate in (H)): 



Po 
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where again uq ~ ^qcd and m denotes the effective heavy quark mass, which differs from the 
perturbative pole mass due to the contribution of gluon condensate. In the nonrelativistic 
rest frame we have^ 
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Then, we can derive that 
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12 \ TT AE^ ^ ' 

Eq.(jH) shows that at (^G^,^) ~ ^qcd ^^e contribution of gluon condensate to the heavy 
quark mass is about Aqcd, i-e. it is linear in the infrared scale of energy, when the operator 
determining this term is of the fourth power in the scale. 

Note, that the second term independent of in the gluon condensate contribution shown 
in results in the correction to the static energy of heavy quark, so that 
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Furthermore, the gluon condensate contributes to ujq in two ways: the first one is explicitly 
given by (0), the second is related with the redefinition of heavy quark mass (m in). 
Indeed, in this case we have to redefine the "large" momentum of heavy quark by the 
substitution for mv by ihv and so on, which means that the resulting change of static 
energy is given by 



Atuo 



m 



m + 5ujn ~ A 
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Then, we can see that after the account for the gluon condensate the displacement of static 
energy can be basically adopted in the mass m. 

Furthermore, we can write down the following relations for the perturbative dependence 
of heavy quark quantities on the scale A: 
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In NRQCD, where |p|/m < 1, the gluon condensate correction to the heavy quark action F tends to zero 
at large virtualities Q — AE as 0{1/Q^) and 0{1/Q^) for the static and dynamic terms, respectively. 
However, the correction remains small even at lower scales. 
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where in the second equahty we neglect the dynamical term and remain the static energy. 
Then the linear dependence on A in m appears in to ways: the first is the direct calculation 
of self-energy diagram for the heavy quark, which results in 

— — = Cmas[A), 
aX 

and the second is contributing from the gluon condensate term due to the AE" dependence 
according to (@)) and (jH)) (the vacuum condensate of gluon operator has the higher power: 
A^), so that 



dX 4 \ TT 

Then, we see that at Ai? ujq the heavy quark mass can be physically independent on the 

dm dm(^) i dm^^ 

dX ~ dX dX 



introduction of factorization scale A, i.e. ^ = ^ + ^ = 0, if 



4_ /« 2 

At i^Gl^) ^ (0.37 GeV)^ @ the evaluation gives 

uo ^ 0.46 GeV. 

Neglecting the dynamical term in the heavy quark virtuality we obtain the following estimate 
of displacement for the heavy quark mass due to the gluon condensate: 

Am ^ ^tuo ^ 0.15 GeV, (7) 

which can serve as the constrain of maximal value, since we expose the minimal virtuality. 

Thus, the main statement on the nonperturbative displacement of heavy quark masses 
remains the following: it is about the confinement scale. However, we can get some definite 
estimates for these values. 

The validity of above consideration is determined by the condition for the formation of 
hadron containing the heavy quark. Indeed, the time for the binding of the heavy quark, 
i.e. for the formation of its wavefunction, in general depends on the hadron contents. So, in 
the heavy-light hadron Hq with a single heavy quark, the static energy for light degrees of 
freedom is of the order of Aqcd, and we get the estimate 

rlHn] ^ 
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which is comparable with the characteristic time for the heavy quark interaction with the 
gluon condensate. In the doubly heavy hadron Hqq, the formation of wavefunction is 
determined by the average size of the doubly heavy system divided by the heavy quark 
velocity 

n-^QQ 2"' 

and it depends on the inverse kinetic energy in the doubly heavy subsystem. So, the calcu- 
lated contribution by the gluon condensate in the heavy quark mass would be inapphcable 
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for the quarks heavier than 20 GeV. However, in the systems composed by charmed and 
beauty quarks the kinetic energy is about Aqcd, and in the reahty we deal with the situation, 
when the effects connected with the formation of wavefunction for the hadron containing the 
heavy quark and the gluon condensate term are competitive. So, for instance, the energy 
shift due to the interaction of coulomb doubly heavy system with the gluon condensate |2| 
is determined by the following expression: 

where is a rational factor depending on the principal and radial quantum numbers n 
and /. Comparing (jSJ with (jl)), we see that despite different approaches these two equations 
can be in agreement with each other, if we substitute of AE ~ En, i.e. if the virtuality is 
determined by the bound energy of the heavy quark in the heavy quarkonium system. This 
fact implies that the heavy quarkonium represents a specific case of the general consideration 
applied to the coulomb system, wherein the virtuality is prescribed to a concrete value. As 
for the numerical estimates, we have to take into account, that one should use (jH) instead 
of (jHl), if the virtuality of the heavy quark in the quarkonium is less that the value following 
from the general form of dispersion law for the quark, i.e. it is less than 0.46 GeV. Otherwise, 
the quark is heavy enough in order to use the coulomb approximation of (jH)). 

To conclude, we have shown that the Operator Product Expansion including the gluon 
condensate results in the following dispersion law for the heavy quark: 

p2 

Po(p) = ^0 + — , 
2m 

where the correction to the heavy quark mass is given by 

A /"^^2 \ 1 

Am = — ( — G,,„ ) 

and the infrared ambiguity in the mass caused by the corresponding renormalon, can be 
cancelled at 

4_ /tts 2 

Of course, the conclusion is drawn to the given, linear order in a^, and the well known 
divergency of heavy quark pole mass with the increase of Og-order probably can be removed, 
if the higher order corrections to the Wilson coefficient of gluon condensate as well as the 
higher condensates will be included into the consideration in the same manner. 
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